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1.  Introduction 


In  the  last  few  years  considerable  theoretical  work  has  been  done 
by  the  groups  at  Flow  Research  Inc.  and  Physical  Dynamics  Inc.  on  the 
effect  of  surface  roughness  on  boundary  layer  transition.  The  basis 
of  this  trork  Is  a phenomenological  theory  of  distributed  roughness 
advanced  by  Merkle,  Kubota  and  Ko  [1].  In  this  theory  the  distortion 
effect  on  the  mean  velocity  profile  caused  by  the  roughness  Is  treated  by  an 
eddy  dlffuslvlty  concept  utilizing  three  empirical  constants.  The  theory 
Is  In  the  same  spirit  as  the  familiar  zeroth-order  closure  methods  (eddy 
viscosity)  of  treating  turbulent  boundary  layer  profllea.  A consequence 
of  this  phenomenological  treatment  Is  that  the  roughness  Is  characterized 
by  a single  parameter,  the  mean  roughness  height. 

Recently  this  phenomenological  theory  has  been  applied  by  Kosecoff 
et  al.  [2]  to  a heated  flat  plate  boundary  layer  In  water  to  determine 
the  effects  of  roughness  and  wall  heating  on  the  neutral  stability  curve. 

One  of  the  conclusions  reached  Is  that  only  a modest  deviation  from  the 
smooth  wall  (Blaslus)  profile  Is  needed  to  drastically  change  the  stability 
characteristics.  Their  results  show  that  the  eddy  dlffuslvlty  roughness 
model  produces  quite  large  deviations  from  the  Blaslus  profile  when 
the  mean  roughness  height  Is  equal  to  the  momentum  thickness. 

The  ala  of  the  present  Investigation  Is  to  determine  whether  the 
magnitude  of  the  distortion  effects  on  the  boundary  layer  profile  predicted 
by  the  phenomenological  theory  of  Merkle  et  al.  Is  reasonable.  To 
•ccoapllah  this  objective  we  will  compare  the  results  of  Kosecoff  et  al. 
with  the  rational  theory  of  Singh  and  Lualey  [3]  applied  to  a flat  plate 
boundary  layer  with  small  amplitude  randomly  distributed  surface  wavlness. 
The  Singh-Uadey  theory  is  a perturbation  approach  for  a quite  general 
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roughness  distribution  which  proceeds  directly  from  the  equations  of 
motion  and  consequently  does  not  suffer  from  empirical  constants.  Since 
it  is  a perturbation  theory,  the  results  are  restricted  to  values  of 
the  expansion  parameter  much  less  than  one.  In  the  present  application, 
however,  this  does  not  seem  to  be  a drawback.  In  contrast  to  the  single 
parameter  characterization  of  roughness  by  the  phenomenological  approach, 
the  Singh-Lumley  theory  shows  that  the  solution  depends  not  only  upon 
the  mean  roughness  height  but  on  the  Taylor  Microscale  of  the  roughness 
distribution  as  well. 

2.  Flat  Plate  Boundary  Layer  with  Small  Amplitude  Waviness 

The  theory  of  Singh  and  Lumley  treats  randomly  distributed  roughness 
described  by  y*h(x,z)  using  a perturbation  expansion  in  terms  of  a 
parameter  e/X  where  E is  the  mean  spatial  roughness  height,  and 

X is  the  Taylor  microscale  of  the  roughness  distribution.  The  assumption 
is  also  made  that  h is  a statistically  homogeneous  isotropic  function  of 
the  spatial  variables.  Characteristic  of  a perturbation  expansion,  the 
expansion  parameter  is  required  to  be  small  compared  to  unity,  i.e., 

f « 1 . (1) 

Physically,  this  is  a rsstrlction  to  roughness  of  small  slope.  The  theory 
then  calculates  by  means  of  matched  asymptotic  expansions,  through 
0(e^/X^),  Che  effect  of  the  roughness  on  the  mean  velocity  profile 
To  obtain  an  order  of  magnitude  estimate  of  the  effects  to  be 
expected  from  the  Singh  and  Lumlsy  theory,  we  consider  a roughness  made 
up  of  a single  wavs  number  a^.  A single  wave  number  simplifies  the 
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calculations  considerably  and  since  the  effects  of  distributed  roughness 
are  expected  to  be  of  the  same  order  as  for  the  olngle  wave  number  case, 
the  added  labor  of  determining  the  spectral  density  for  distributed 
roughness  does  not  seem  worthwhile.  A roughness  characterized  by  a single 
wave  number  corresponds  to  small  amplitude  wavlness  distributed 
randomly  and  Isotropically.  This  type  of  wavlness  Is  commonly  encoiintered 
in  practice  and  results  from  the  manufacturing  process. 

The  mean  velocity  in  the  boundary  layer  In  the  x (streamwlse)  direction 
is  given  by  Singh  and  Lumley  as 

<p\yiy)  “ , (2) 


where  u(y)  Is  the  velocity  distribution  for  the  smooth  wall  case  and 


Is  the  perturbation  due  to  surface  roughness  and  Is  found  to 
av 

be  given  by 


(R  Up  5)^^^  E(a)  (Ks)da 


where 


(3) 


'•'->1/3 
a - y(R  UqO) 


f 


r 


4(s)  ■>  influence  function  tabulated  by  Singh  and 
Lumley  in  their  Table  I, 

Uq  • slope  of  the  undisturbed  velocity  profile  at 
the  well,  y"0,  and 

1(5)  ■ roughness  spectral  dsnslty. 
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The  roughness  spectral  density  E(a)  Is  defined  by 


E(a)da 


For  a roughness  characterized  by  a single  wave  number,  the  spectral 


density  Is  given  by 


where  5(a-a^)  Is  the  Dirac  delta  function. 

The  expansion  parameter  e/X  may  be  related  to  the  roughness  wave 
number  as  follows:  By  definition 


In  the  present  case,  the  surface  Is  described  by  y>h(x).  Independent  of 


z,  and  hence 


E(5)da 


For  discrete  wave  numbers  for  which  the  spectral  density  Is  given  by 
Bq.  (5),  Eq.  (6)  yields  the  simple  relation 


The  spectral  density  given  by  Eq.  (S)  also  allows  Bq.  (3)  for  ^2 
be  evaluted  by  quadrature  yielding  the  result 


with  the  eld  of  Eqs.  (5)  and  (7)  the  mean  velocity  given  by  Eq.  (2) 


We  note  that  Eq.  (9)  la  written  In  terms  of  dimensionless  variables  with 


reference  quantities  as  yet  unspecified 


the  basic  flow  Is  chosen  to 


be  the  Blaslus  flat  plate  solution  which  can  be  written  In  terms  of 
similarity  variables  as,  cf.  Rosenhead,  p.  222,  [4], 


where  the  similarity  variable  n !•  related  to  x,  y and  R by 


In  Eq.  (9)  Is  given  by 


For  a Blaslus  profile  Eq.  (9)  becomes 


and  from  Rosenhead,  Table  V.l,  f, 


0.46900 


To  compare  with  Ref.  (2),  we  take  as  the  reference  length  the 


boundary  layer  momentum  thickness  6.  The  Reynolds  number  is  based 


on  6 and  y as  given  by  Eq.  (11)  is  made  dimensionless  by  9.  The 


present  theory  is  therefore  characterized  by  the  three  parameters  R, 


Since  6 is  taken  as  the  characteristic  length,  x in  Eqs.  (11)  - (13) 


is  not  arbitrary  but  is  determined  by  the  definition  of  momentum  thickness 


which  in  Blasius  variables  is 


From  Rosenhead,  Table  V.l 


f'(l-f’)dTi  - 0.46960 


and  hanca  the  relation  between  x/9  and  R.  from  Eq.  (14)  ia 


2.267325  R, 
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* Furthermore,  the  relation  between  y/6  and  n,  from  Eq.  (11),  now 
reads 

9 

0 " 0.46960  ’ 

3.  Results 

Kosecoff  et  al.  [21  In  their  calculations  of  flat  plate  mean  velocity 
profiles,  the  results  of  which  are  plotted  In  their  Fig.  1,  use  values 
of  R0~175  and  c/6>0.3,  1.0  and  2.0.  The  curve  for  e/6*0.3  Is  Indistinguishable 
from  the  Blaslus  profile,  however,  significant  departures  occur  for  the 
cases  e/8-1  and  2. 

* Calculations  were  performed  for  the  Singh  and  Lumley  theory  using 

Eqs.  (12)  and  (13)  together  with  Table  I of  Ref.  3 for  Rq-175  and  c/0-1 
> ® 

and  2.  Values  of  the  roughness  wave  number  were  chosen  In  the  range 
0.05  to  0.20  keeping  In  mind  the  restriction  ca^«l.  Comparison  of  the 
mean  velocity  profiles  for  the  two  theories  are  shown  In  Figs.  (1)  and 
(2)  for  e/0-1  and  2 respectively.  To  show  more  clearly  the  differences 
In  predictions  of  the  two  theories,  the  disturbance  velocity  ratio 
6u/u  Is  presented  In  Figs.  (3)  and  (4)  for  the  same  values  of  e/0.  The 
latter  two  figures  show  quite  clearly  that  the  Singh  and  Lumley  theory 
predicts  a velocity  profile  disturbance  nearly  an  order  of  magnitude 
smaller  than  the  phenomenological  theory  of  Merkle  et  al. 

Wa  note  that  In  Figs.  (2)  and  (4)  the  curves  for  a^-0.1  are  beyond 

* the  limits  of  vslldlty  of  the  present  theory  because  ea^  is  no  longer 
much  less  than  one.  Also,  the  nagstlve  velocity  predicted  by  the  Singh 
and  Lumley  theory  near  the  wall  Is  evidence  of  the  theory  breaking 
down  In  this  region  becauae  of  llnaarlslng  the  boundary  condition  there. 

\ 
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The  reason  for  the  large  difference  In  the  predictions  of  the 
Kosekoff  et  al.  theory  and  that  of  Singh  and  Lumley  is  almost  surely 
not  the  fact  that  one  is  ad  hoc  while  the  other  is  exact,  but  rather 
the  fact  that  the  former  is  calibrated  against  typical  sand  roughness 
drag  coefficients,  corresponding  to  e/X  » 1,  while  the  latter  is  valid 
for  c/X  « 1.  Hence,  the  theories  are  not,  properly  speaking,  comparable. 
The  theory  of  Kosekoff  et  al.,  however,  particularly  as  regards  its 
stability  predictions  in  the  presence  of  heating,  is  being  put  forward 
as  relevant  to  small  axlsymmetrlc  underwater  vehicles.  The  type  of  rough- 
ness present  on  such  vehicles  is  due  to  the  manufacturing  processes  used, 
and  since  every  effort  is  made  to  produce  a smooth  finish,  it  is  most 
likely  to  be  a small-slope  wavlness,  rather  than  a large-slope  sand 
roughness.  Hence,  we  would  expect  the  theory  of  Singh  and  Lumley  to  be 
applicable  to  these  bodies,  and  consequently  that  a given  roughness  height 
would  be  far  less  detrimental  than  the  theory  of  Kosekoff  et  al.  would 


suggest. 
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